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Abstract. — We consider the Klein-Gordon equation (KG) on a Riemannian surface M 

dtu - Am - m 2 u + u 2p+1 = 0, p G N*, (t,i)6RxM, 

which is globally well-posed in the energy space. This equation has an homoclinic orbit to the origin, and 
i i in this paper we study the dynamics close to it. Firstly, using a strategy from Groves-Schneider, we show 
that there are many solutions which stay close to this homocline during all times. We point out that the 
solutions we construct are not small. 
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^ 1. Introduction, statement of the main results 

1.1. General introduction. — Denote by M a compact Riemannian manifold without boundary 
of dimension 1,2 or 3 and denote by A = Am the Laplace-Beltrami operator on M. In this paper we 
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are concerned with the following nonlinear Klein-Gordon (KG) equation 

J <9 2 u - An - m 2 u + u 2p+l = 0, (i, x) G E x M, 
| u(0,x) = -uo(x), «9 t n(0,a;) = ni(x), 

where p > 1 is an integer, and (u ,ui) G H X (M) x L 2 (M) are real-valued. 



(KG) 



It is well-known that there exists a Hilbert basis of L 2 (M) composed with eigenfunctions (e n ) n >o 
of A. Moreover, there exists a sequence = Xq < X± < ■ ■ ■ < X n < . . . so that 

-Ae n = X 2 n e n , n > 0. 

In the sequel, we define the scalar product on L 2 (M) by (/, g) = — — j fg, where \M\ denotes the 
volume of M, we assume that ||e n ||^2 = 1 and we set cq = 1. 

We make the assumption 



Assumption 1. — In (KG), m satisfies 

< m < Ai. 

In the particular case M = S 1 we prefer to index the basis on Z and we set 
(1.1) eo(x) = 1 e n (x) = V2cos(nx) and e_ n (x) = v / 2sin(na;) for n > 1, 

and we have Ai = 1. 



We make the following assumption on M and p G N* in ( KG ) 
Assumption 2. — The manifold M and the integer p satisfy either: 

• M = S 1 = M/2ttZ andp>\ 

• M is any compact manifold without boundary of dimension 2 and p > 1 

• M is any compact manifold without boundary of dimension 3 and p = 1. 

The stationary solutions of (KG) (solutions which only de pen d on the space variable) are exactly 
the constants u = 0, u = m 1 ^ and u = — m 1//p (see Lemma 2.1). The origin is an equilibrium with 



an unstable direction. In fact, the eigenvalues of —A — m 2 are the (A| — m 2 )k<^n- Since < m < X\, 
the case k = only, gives the hyperbolic directions, corresponding to the solution exp(mt) for t > 



this work we aim to study the dynamics of (KG) near this particular trajectory. 



(resp. exp(— mt) for t < 0). It turns out that (KG) admits an homoclinic orbit to the origin and in 



In the case M = S 1 we can precise the dynamics around the equilibrium u = m 1//p (the study 
near u = —m l / p is similar since the non-linearity is odd). We linearize the equation near this point 
(u = m l /' p + v) and we are lead to study the spectrum of (— d 2 + 2pm 2 )v. This operator on S 1 is 
self adjoint and has pure point spectrum j 2 + 2pm, j G N. Using the Birkhoff normal form theory, 
Bambusi [1] (see also |3]) has shown that, for a generic choice of m (in order to avoid resonances 
between the frequencies y j 2 + 2pm, j G N), the solutions of (KG) with initial datum of the form 
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m Vp _|_ Vq with vo small enough, remain close to the equilibrium point m 1//p for very long time. See 
also the work |8j of J.-M. Delort for quasi-linear equations. 

Moreover, using KAM theory, C. E. Wayne |26) and J. Poschel |24| have proved that, for a generic 
choice of m, there exist many quasi-periodic solutions near w = m 1//p . See also the recent work of 
Berti-Biasco-Procesi [4] for derivative wave equations. Observe that quasi-periodic solutions can be 
constructed even if the equilibrium has a finite number of hyperbolic directions (see [5]). 

In higher dimension and for a general manifold, few is known. In the case M is a Zoll manifold, 



Bambusi-Delort-Grebert-Sjeftel [2] have developed a Birkhoff normal form theory for (KG) near an 
elliptic equilibrium. Up to now, there is no KAM-type result for (KG) in dimension greater than two. 
In this work, we describe some possible behaviours of (KG) near the homoclinic orbit. We state 

rigin to turn around t 
and its Corollary 1.3). We stress 



the existence of solutions that travel from a neighbourhood of t he o rigin to turn aro und the equilib 



1.2 



num m" /f/ , close to the homoclinic connection to (see Theorem 
out that these results do not require non resonance conditions. In Theorem |1 .5 [ we state the existence 
of some large solutions connecting periodically a neighborhood of size r\ to itself, doing some loops 
close to the homoclinic orbit during a time of order (ln(^)) 2 ~ <5 . 

The existence of homoclinic or heteroclinic connections to periodic or quasi-periodic solutions is a 
question of interest : for example, in the case of parabolic PDE's on a one dimensional bounded domain, 
homoclinic orbits to equilibrium points can exist but there cannot exist any homoclinic connections to 
periodic orbits (see [17J). In some other cases some homoclinic connections to small periodic orbits 
exist while homoclinic orbits to do not exist or remain an open question (see for instance |18[ I19|. [16J). 
More generally, such homoclinic connections to the center manifold often appear when an homoclinic 
orbit does not persist after a perturbation of the system (see for instance |25j ). We refer also to [6j, 
where the existence of homoclinic solutions to a beam equation is studied. 

In this paper, we get the existence of a large family of homoclinic connections to the center manifold 
(more precisely heteroclinic connections to some solutions lying in the center manifold). The proof 
is based on a perturbative method which is classical (see for instance |15|, I18| ) for finite dimensional 
reversible systems (systems which anticommute with a symmetry S) : the key idea is that if a solution 
hits the reversibility plane {u \ u = Su}, then the latter solution is reversible. This method requires 
more computations in the infinite dimensional cases. It was already performed in a situation close to 
ours by M. Groves and G. Schneider 1121 113] : but in their case, they work in the neighbourhood 
of a bifurcation and get small homoclinic solutions, while in our case the size of the solutions is of 
order 1 and this requires some additional work on the linearised system. 

Finally we also mention the recent book of Nakanishi & Schlag [22j on invariant manifolds in the 
context of dispersive Hamiltonian PDEs. 



1.2. A first motivation for studying (KG), 
equation reads 

(1.2) d?w-Aw + w-f(w) 



Recall that the usual non linear Klein-Gordon 



0, (t, x) € K x M, 

where / is a non linear function. Let us show that if there exists a nonzero equilibrium, the equation 



near the smallest equilibrium is of the form ( KG ) (but with a general non linearity) . The equilibrium 
of (1.2) are the real constants wq satisfying /(u>o) — wq = 0. Observe that if one performs the change 
w — Wo, then there exists a non linear function g such that u satisfies 

3 2 t u -Au- (f(w ) - l)u - g(u) = 0. 



of coordinates u 
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Given that / is non linear, /'(0) = and /'(0) — 1 < 0. Thus, if / is C 1 and if there exists an 
equilibrium wq ^ 0, then by the intermediate value theorem, we get that the smallest equilibrium wq 



satisfies f'(wo) — 1 > 0. This means that near the smallest non zero equilibrium of (1.2), the equation 
is of the form 

d 2 u — An — m 2 u — g(u) = 0. 
Here we study the particular case of a non linearity g of the form u 2p+1 . This case arrises for m 2 



(which belongs to (0, 1)), wq 



-m» and f(w) = (w + m? ) p — m p — (2p + l)m w. 



2p 



1.3. Another motivation of the problem. — Here we are inspired from the works [11, 112L 113] 
of M. Groves and G. Schneider. Consider the non linear wave equation on the real line 

(1.3) d?w -d 2 x w + w- w 2p+1 = 0, (t, x) G R x R. 



One can be interested in the construction of "modulated pulse solutions " for (1.3) which are solutions 
of the form 



(1.4) 



,x — fit t — [3x^ 
w(t,x) = u{— — , j. 



Cj3 C/3 

where u(s,y) is a 2ir— periodic function in y, and where (3 € (0, 1) and eg > 0. In the particular case 
cp = (1 — /3 2 ) 1//2 , (1.4) is the Lorentz transform which preserves (1.3) and in general we get 

„2 



;i.5) 



d 2 u - d 2 u 



1-p 



(u-u 2p+1 )=0, (s,y)£RxS 1 . 



Now let m E (0,1) and /3 G (0,1), then we can choose cp > so that ^pi = fn 2 £ (0,1), which 
is (KG) up to the change of unknown u i — )• u l l p w. As a consequence, for each m G (0, 1), any solution 



to (1.5) provides a one-parameter family of solutions to (1.3). 



1.4. Hamiltonian structure of (KG) 



(1.6) 



H 



V x u\ 2 + v 2 



Denote by v = dfU and introduce 
1 



2 1 

m u I + 



2p + 2 



,2p+2 



Then, the system (KG) is equivalent to 
(1.7) 



6H_ 

8v 



5H_ 

5u 



We write 



where 



i(t, x) = a n (t)e n (x), v(t, x) = ) j b n (t)e n (x), 



n=0 



n=0 



(a n ) n( z N G ^(N,^) := {x = (x n )„ eN | ||a;||^i = ^(1 + \l)\x n \ 2 < oo}. 
(ftn)neN G £ 2 (N,R) := {x = (x n ) neN | ||x|| 2 2 = ^ |x n | 2 < oo}, 



neN 
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in such a way that to the continuous phase space X := H 1 x L 2 corresponds the discrete one W h 1 x I 2 . 
We endow this space with the natural norm 

ll^ll = INI./? 1 + II^IIl 2 ) f° r x = (u,v). 

We also define the two dimensional subspace X^ C X spanned by the stable/unstable modes (hyperbolic 
modes) of X 

(1.8) X h = {X € X | Vn > 1, a n = b n = 0}, 

as well as the infinite dimensional subspace X c C X spanned by the center modes (elliptic modes) of X 

(1.9) X c = {X e X | a = b = 0}. 

We clearly have X = X^ © <Y C . Denote by P : <Y — X^ the projection onto Af^ and define Q = 
I — P : X — > X c the projection onto X c . In the sequel, for X £ X, we use the notation X^ = PX and 
X c = QX. 



In the coordinates (a n ,b n ) n >o, the Hamiltonian in (1.6) reads 



(1.10) 



H 



i 00 



m )a n + b n 



n=0 



+ 



1 

2p + 2 



M 



2p+2 



k=0 



and the system (1.7) becomes 

' a n = b n , n > 

(1.11) 



6„ 



-(X n -m )a r , 



r / ~ n2 P +i 
/ > afcefc(a;) e n (x)dx, n > 0. 



:i.i2 



1.5. The homoclinic orbit. — The space-stationary solutions of (KG) exactly correspond to the 
solutions of ( |1.11[ ) satisfying a n = b n = for n > 1. In this case, the equation on (ao, 6o) reads 

' a = b 

U 2 2p+l 

t)o = m oo - a , 

and this system possesses an homoclinic solution to that we will denote in the sequel by h : t \— > 
(a(t), j3(t)) and we will denote by /Co the curve which is described (see Figure [T]). Indeed we can 
explicitly compute 

smh.(pmt) 



m 1 /p(n + 1 W(2p) 



-m 



l/p+i 



(P+1 



,l/(2p). 



COS 



( cosh(pmt)) 
and we have the bounds 

(1.14) |a(t)| < Ce- m l*l, |/3(t)| < Ce~ m ^, for all tel. 



h(pmt)) 



1/p+i ' 



For 77 > denote by /C^ the trajectory of (1.12) given by the initial conditions ao(0) = 77, 6q(0) = 

1). This orbit is periodic with a period of order In—: most of that time is dedicated to 

U *1 

cover the very principle and the very end of the loop, i.e. when (ao, bo) is close to (77, 0). 



(see Figure 



1. In the particular case M = S 1 , according to definition the index n cover Z in the previous sums. 



BENOIT GREBERT, TIPHAINE JEZEQUEL AND LAURENT THOMANN 




FIGURE 1. Phase portrait for the space-stationary set a n = b n = for n > 1. 
1.6. The main results on the Klein-Gordon equation. — Under Assumptions [T] and [2j our 



first result is that the equation (KG) is globally well-posed in X = H (M) x L 2 (M). 



Theorem 1.1. — Let M, p and m satisfy Assumptions^ and^ Then for all (uo,Wi) £ X there exists 
a unique solution to (KG) 

u £ C°(K; iJ-^M)) DC 1 (R;L 2 (M)). 
Moreover, u is bounded in the energy space: sup ||(u, < C. 



The proof of Theorem |l.l| is very simple: it relies on a fixed point argument in H 1 (M) combined 
with the conservation of the energy. 

Our next result describes possible dynamics near the homoclinic orbit: We show that there are many 
solutions which stay close to h. The proof of this result is inspired from the work [llj of M. Groves 
and G. Schneider. In our case, the novelty is that we have not to assume m to be small, i.e. we are 



not dealing with small solutions of ( KG ) . 



For any function / £ L 2 (M), denote by / c = / — / /. Then we can state 

Jm 

Theorem 1.2. — Let M, p and m satisfy Assumptions^ and^ There exist C > and e$ > so 
that the following holds: Let (f,g) £ H (M) x L 2 (M) be such that WfcWn 1 + Hs'cIIl 2 — £ 0; then there 



exists a solution u to (KG) so that 

(1-15) d t u(Q,-)=g c , n c (0,-) = / c , 

and so that for all t > 

(1-16) \\u(t,-)-a(t)\\ H1 + \\d t u(t,-)-p(t)\\ L2 < C{\\f c \\ m + \\gcWv). 



Furthermore, as t — > +oo, u(t, •) tends to a solution of the local center manifold W° (see Defini- 
tion 4-6). 
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Notice that Theorem 1.2 is not a stability result since we only fix the elliptic part of the initial 
datum, and the statement of the theorem is that there exists a way to define the hyperbolic part such 
that the solution satsifies (1.16). For example if f c = g c = 0, the corresponding u is the homocline h. 

The equation ( KG[ ) is reversible. More precisely, the equation (KG) is invariant under the transfor- 
mation u(t,x) i — > u(—t,x). This corresponds to the symmetry S defined by 

S(a n ) = a n S(b n ) = -b n , for n > 0. 

This is a symmetry of reversibility: Namely, denote by Vh the vector field (1.11 ), then Vh°S = — SoV H . 



Using this symmetry we are able to precise the result of Theorem 1.2 in the case g = 0. 

Corollary 1.3. — Let M, p and m satisfy Assumptions^ and [£| There exist C > and £q > 
so that the following holds: Let f £ ff 1 (iVf) be such that ||/ c ||#i < £o> then there exists a solution u 
to fKG[ ) so that 

(1.17) d t u(0,-) = Q, «c(0,-) = / Cl 

for allteR 

u{-t) = Suit), 

\\u(t, •) - a (t)\\ H i + \\d t u(t, •) - (3(t)\\ L 2 < C\\f c \\ m . 
+oo, u(t, •) tends to a solution of the local center manifold W c . 

> +oo, then u tends to Sw when 



and so that 
Furthermore, as t 



Since u is reversible, observe that if u tends to w when t 
t — > — oo. 



Remark 1.4- — Corollary 1.3 proves the existence of heteroclinic orbits from some orbit on the center 
manifold W c to another one. It would be interesting to precise the dynamics on the center manifold. 
Actually we can expect the existence of invariant tori on W c which are perturbations of the tori 
obtained for the linearized equation around the homoclinic orbit (see section 3). Nevertheless, to prove 
this, we would need a KAM theory for infinite dimensional tori which is an open problem up to now 
(see however [7j)- Notice that even if such theory would exist, it would prove that most of the tori, 
but not all of them, of the linearized system are preserved. Thus it would be impossible to conclude 
that our heteroclinic orbits tends to KAM tori of W c when t — > ±oo as we can expect. 



Our next result is a description of another possible dynamics in (1.11) (and thus of (KG)): There 
exist trajectories which are close for a long time to /C^, turning around the equilibrium point. Indeed we 
prove that if the initial datum is close to the origin in the mode (precisely (ao = rj,bo = 0) with r/ > 
small enough) and even closer to the origin in the mode n, n > 1 (precisely a n ,b n = 0{t] 2 )) then the 



solution will describe in the mode several loops close to the homoclinic orbit of (1.12) during a time 
of order — In 77 while the mode n will describe even more loops around the origin. Roughly speaking, 



this solution follows approximatively the periodic orbits of the two-dimensional system (1.12) 



Theorem 1.5. — Let M , p and m fixed and satisfying Assumptions^ and^ There exists r/Q > such 
that if < 7] < t]q and if the initial datum satisfies ao(0) = T), bo(0) = and X^n>i(-^n a n + ^n) — V 
then there exists C > so that 
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(i) There exists T v > such that -q In - < T v < Cln -, which satisfies 

bo(T v ) = and \ao(T v ) — rj\ < Crf . 

(ii) Let 5 > 0. The solution of ( l.llj ) satisfies for all \t\ < ( In ^) 2 5 

^((a (i)A(i)),*V) < Cr? 5 / 2 and £ (A 2 a 2 (*) + b 2 n (t)) < Cifl 2 . 



n>l 



Remark 1.6. — The time T v is in fact the time necessary for (ao,bo) to complete one loop around 
(m 2 ,0), and the previous result shows that the solution remains rf— close to ]C ri during 
loops. 



In 



1.7. Plan of the paper. — In Section [2] we prove Theorem LI 



of the linearisation of (KG) around the homoclinic orbit (Theorem 



Section [3] is devoted to the study 
O. In Section [| we construct the 



center manifold for the Klein-Gordon equation and prove Theorem 1.2 In the appendix, we recall a 



result which is useful in the proof of Theorem 3.1 



Notations. — In this paper c, C > denote constants the value of which may change from line to 

line. These constants will always be universal, or depend on the fixed quantities m and p. 

We denote by N the set of the non negative integers, and N* = N\{0}. We set X = i^ 1 (M) x L 2 (M). 

Acknowledgements. — The authors want to thank Eric Lombardi and Romain Joly for clarifications 
in dynamical systems. T.J. also thanks Guido Schneider for his kind invitation to Stuttgart. 



2. Preliminaries: Proof of the first results on (KG) 



2.1. Stationary solutions of (KG). — We look for the time-independent solutions to (KG), which 
correspond to equilibrium points for the system ( |1.11[ ). 



Lemma 2.1. — The stationary solutions to ( |KG[ ) are u = 0, u = m l / p and u = —m l l p . 

Proof. — Assume that u(t, x) = g{x) S ^(M) is solution to ( |KG[ ). We apply V x to the equation and 
get 

-AVg - m 2 Vg + (2p + l)g 2p Vg = 0. 
We multiply this equation with Vg and integrate (by parts) on M 

(2.1) / {Ag) 2 + (2p + l) [ g 2 P\Vg\ 2 = m 2 [ \Vg\ 2 . 

JM JM JM 

oo 

Now, write g = '^^gk&k- Since (\k)k>o is non decreasing and < m < X± we have 



fc=o 



„r i \Vg\' 

• M 



m 



oo oo „ 

J2xl 9 l<J2Mgl= / (Agf 
fc=i fe=i Jm 



This inequality together with (2.1) implies that Vg = 0, and we conclude. 



□ 
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2.2. Proof of Theorem 1.1 

(2.2) 

oo 

and for u = a p e p we define 
p=0 



To begin with, we set 



J 



^ oo 

fc=l 



2 2\ 2 , 1.2 

k -m )a k + 6 fc 



(2.3) 



U := u r 



(1 - n )u = y^a P e p 
P =i 

the spectral projection away the mode 0. In the sequel we will need 
Lemma 2.2. — Let u £ H l (M). Then 

• When M has dimension 1 or 2, for all 2 < q < +oo 

||f ||i?(Atf) < CgJ 1 / 2 . 

• When M has dimension 3, for all 2 < q < 6 

||tf|U?(M) < CqJ 1 ! 2 . 

Proof — By Sobolev, in each of the previous cases, there exists C q > so that for all U 6 H l (M) we 
have ||J7||i«(Af) < C<jll^ r ll_H" 1 (Af) an d the result follows. □ 

We then can prove the following a priori estimate 



Lemma 2.3. — Let u(t,-) be a solution of (KG) and denote by H° the constant value of H(u(t, ■)) 
along this trajectory. For all t for which the solution is defined, we have 



,. oo 

(2-4) / {{d t uf + \Vu\ 2 ) = bl + £ (bl + Xlal) < 2H° + — ^— - 

Furthermore there exists C = C(H°) > such that |ao(i)| < C for all t. 



m 



2+2/p 



Proof. — We assume that M has dimension 1 or 2. By the conservation of the energy for (KG), we 
have 



/ 

JM 



, ((3 t u) 2 + |Vu| 2 ) + -l T 
We apply the Young inequality 



u: 



2p+2 



2H° + m 2 



M 



u 2 . 



Al 



(2.5) 



cic 2 = (eci)(e 1 c 2 ) < —c\ + -^—c r 2) ci, c 2 > 0, e > 0, - + - = 1, 



with ci = u 2 , C2 = 1, ^ = p + 1 and e = m 2 /(p +1 ) and deduce (2.4). 
Recall the notations (2.2) and (2.3). Then we have 

(2.6) 



By (2.4), &o an d J are bounded and moreover by Lemma 2.2 and ( |2.4[ ) we deduce that J M |£7| 9 d2; < 
CJi/ 2 for all q > 2. Thus (|2~6|) gives H° as a polynomial in qq with bounded coefficients and thus oq 
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has to be bounded by a constant depending only on H°. 

The argument is similar when M has dimension 3 and p = 1. □ 

Remark 2-4- — It is easy to check that the points u = ±m 1//p correspond to the minimum of H, and 

that for this choice H° = — - — -m 2+2 l p . In particular, the r.h.s of (2.4) is nonnegative. 

2(p+l) 1 — r 



Proof of Theorem 1.1 — First we show that (KG ) is locally well-posed. Let (uq, U\) 6 X and for T > 
define the space 

E T = {ueC°([-T,T};H 1 (M))nC 1 [[-T,T};L 2 (M)), sup \\(u,d t u)\\x < 2\\(u , Ul )\\ x }. 

\t\<T 

Denote by K(t) the free wave propagator defined by 

/ \ / \ / / — r\ sux(t\f— A) 
K (t)(uo, Ul) = COS (W— AlMn H ; — ui. 

Then we show that the application $ defined by 

*(t)u = iT(t)(n , ui) + /' Sln(t "j/^ ( m 2 n _ u 2 P+ i) (s)dfi) 

Jo v— A 

is a contraction of for T small enough. Namely, by the Sobolev embeddings and the fact that 
1 1 -^11 A") — we obtain that we can take T = comin (m -2 , ||(^0; u i)||^' 2p ) where cq > is a small 
absolute constant. 

In order to prove the global well-posedness, we iterate the previous argument. We then obtain a 
sequence of times T n with T n +i = Co min (m -2 , \\(u(T n ),dtu(T n )\\~^ 2p ^ such that the solution is de- 
fined on (- Y^7=o T ™> X^°=o T ™)- B Y Lemma g for any t > 0, \\(u(t), d t u{t))\\ x 2p > C(H°), hence 
Sn>o Tn = +oo which gives the result. □ 

Remark 2.5. — The result of Theorem 1 1 . 1 1 indeed holds for any m£R. 



3. The linearized equation around the homoclinic orbit 

We first consider the simplified model where we omit every interaction term containing a product 
of two modes a^, ai with k,£> 1. The idea is to capture the principal dynamics when we are close to 
a n = b n = for n > 1. The resulting toy model reads 

b 



(3.1) 
and for n > 1 
(3.2) 



a 



9Q 



m 2 oo 



2p+l 



-(\ z n - m z )a n - (2p + l)a 2 Q P a n . 



Indeed, if one considers that ao is given by (3.1), the system (3.2) is the linearisation of (KG) around 
the solution u = ao- In particular, denoting V(t) = (2p + l)aQ P (t) and w = Ylk>i a k e k, (3.2) is 
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equivalent to the following linear wave equation with time-dependent potential 
(3.3) d?w -Aw- m 2 w + V(t)w = 0. 



Notice that the system (3.1), (3.2) is not Hamiltonian for the canonical structure, but it is reversible. 
On the other hand, for ao given, the system (3.2) is Hamiltonian for the canonical structure. 
We prove that such an infinite dimensional system has homoclinic orbits to invariant tori of arbitrary 
large dimension (finite or infinite). These tori correspond to the case oo = &o = which is stable 
by (3.1). In that case (3.2) becomes the standard harmonic oscillator in infinite dimension with 
frequency vector oj(m) = (oj n ) n >i and uj n = \/X^ — m 2 . Namely, the tori are given by 

% := {{a h ,b k )h>\ | (X 2 n -m 2 )a 2 n + b 2 n = c 2 n }, ce£ 2 (N,IR), 

and we have 

Theorem 3.1. — Let (ao(0), ?>o(0)) £ an d (a n (0), b n (0)) n >% G h 1 x £ 2 . Then the corresponding 
solution of (3.1) 7 (3.2) is homoclinic to T c for some c £ £ 2 (N, M). 

This result can also be interpreted as a linear scattering result for the system (3.2). 

Proof. — Take the homoclinic (a, j3) in the plane (ao,bo), and observe that / a 2p (t)dt < oo. Now, 

Jr 

for n > 1 we introduce the norms on IR 2 , \x, y\ n = ((A 2 - m 2 )x 2 + y 2 ) 1/2 and the norm on h x I 

that the trajectory 



, b n ) n >i || = ( \ a m b n \n) ■ Then for all n > 1, we deduce from Lemma 



n>l 



A.l 



(a n , b n ) satisfies \a n , b n \ n < C\a n (0), b n (0)\ n . Next, we denote by S the flow of the 
i.e. with ao = 0. Then the solution of (|3.2|) reads 



(3.4) 



a n (t) 

b n (t) 



S(t) 



On(0) 
&n(0) 



(2p+l) ts(t- S ) 
Jo 







a 2p (s)a n (s) 



inear part of (3.2 ), 



ds. 



By construction S preserves the norm . | and thus 



+oo 



S(-s) 







a 2p {s)a n (s) 



ds 



a 2p (s)\a n ( S )\ds 



(3.5) < C|a n (0),6 n (0)| n , 

where we used the boundedness of a n and the integrability of a 2p . Denote by 



b+ 



therefore by (3.4) we get 

a n (t) 
bnit) 



b n (0) 



S(t) 



r+oo 

(2p + l) S( 
Jo 





a 2p (s)a n (s) 



ds, 



a,,\ 
bt 



\ r+oo 
) n < (2p+l) J a 2p (s)\a n (s)\ds 



< C\a n (Q),b n (0)\ 7 



+oo 



a 2p {s)ds, 
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which in turn implies 



a n (t) 

bn(t) 



S(t) 



0, when t 



+00. 



lives on the torus T c with c 2 



(A 2 -m 2 )(a+) 2 + 
□ 



Finally we remark that the trajectory t 1— > S(t) 
(6+) 2 for n e N. 

Remark 3.2. — If w(m) is non resonant (for example when M = S , this happens for a generic choice 
of the mass m, see [3]). S(t)(^a^, bn) t>0 densely covers the torus T c . 



Remark 3.3. — Notice that in the proof of Theorem 3.1 we crucially use that ao = a, i.e. that we 



linearize around the homoclinic orbit. A natural question is to ask whether this resul t sti 
(ao, bo) is any periodic solution of (3.1). In that case J K aQ P (t)dt = +00 and Lemma 



A.l 



1 holds if 
does not 

apply. In particular it is not clear at all that the trajectories remain bounded. This seem to be a 
difficult and interesting problem related to the reducibility of (3.3) when V is not a priori small (the 
case V small and quasi-periodic in time can be solved by a KAM approach, see [9j 110] for the last 
results in the Schrodinger case). 



We end this section with two results on the linear evolution of (3.2), which will be useful in the 



proof of Theorem 1.2 Introduce the notations 



Z 



d t z 



then (3.1), (3.2) can also be written 



1 

A + m 2 - (2p+ l)a 2 P(t) 



d t Z = C l Z. 



The next result is a direct consequence of Theorem 3.1 

Proposition 3-4- — The restriction £* : X c — > X c generates a two parameter group K(t,r) so that 



sup \\K(t,T)\\ LM) <C. 

t.r6K 



In other words, every solution to dtZ c 



C C Z C is bounded. 



Denote by C 
the space Xh = 
Lemma 3.5. 



\x h 







1 



m 



(2p+l)a 2 P(t) 



the restriction of £* to Xh- We can identify 



I 2 and we denote by ( , ) the Euclidian scalar product. Then 
The equation 



dtZh = £-hZh 

admits two solutions a(t) and p(t) defined on M + so that 

(3.6) <r(0)=(j), p(0) 
(see Figure^ and 

(3.7) \a(t)\ < Ce 



-rat 



\p(t)\ < Ce m \ for all t > 0. 
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The dual basis {cr*,p*} to {a, p} in satisfies 

(3.8) \o-*(t)\<Ce mt , \p*(t)\ < Ce~ mt , for all t > 0. 



plane {X = SX} 




Figure 2. Linearized system at t — and reversibility plane. 



Proof. — Let p, € K, then 1 1— > p 



a(t) 
a(t) 



/' i p^t) ^ * S a so ^ u ^^ on °^ ®t^h = ^h^h Then ill YH'W 



of (1.13) we deduce that this solution corresponds to a and we obtain the first bound in (3.7). With 



7 



7 



a suitable choice of p we get (3.6). 
To find p, we write p = 



, and look for 7 of the form 7 = az, when t > 1. Then z satisfies the 



equation az + 2cii = and therefore 

/t 
(d(r)) _2 dr, Vt>l. 

It is then straightforward to check that 7 = az can be extended to a C°° function on M. Then 
the bound on a gives \z(t)\ < Ce 2mt and \z(t)\ < Ce 2mt for all t > 0, which in turn implies that 
17(01 < Ce mt and |7(i)| < Ce m * for all t > 1. Hence the second bound in (3.7). We can moreover 
choose the initial value p of the form claimed in the statement of the lemma. 



Then we can explicitly compute 



1 



a 



7 



1 



P 



/?7 - 7/? V -7 J ' ' /3 7 - 7/3 V Z 3 
where ^7 — 7/3 is a non vanishing constant, since it is the Wronskian of the equation 

y - m 2 y + (2p + l)a 2p y = 0. 



The bound (3.8) then follows from (3.7). 



□ 
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4. Homoclinic orbits to the center manifold 



In this section we prove Theorem 1.2 



As we mentioned in the introduction, this part closely follows 
the argument of Groves-Schneider |13] , Nevertheless in contrast to their work, we are not dealing with 
to small and thus with small solutions. Even if some results already appear in |13] . we reproduce here 
all the proofs for the convenience of the reader. First we recall the general strategy, illustrated by 
Figure [3] 

1. To begin with, we consider a spatially truncated system for which the non-linearity sees only the 
elliptic (or central) modes of small amplitude. We construct global solutions to this system that 
are close to the homoclinic orbit, namely the hyperbolic part is close to h for all times t > 0, and 
the elliptic part is small only for t < e _1 where e is a small parameter. (Subsection 4.2) 



In the next steps we will prove that the elliptic part of these solutions actually stay small for all t, 
i.e. that the set of all these solutions is the global center-stable manifold W cs of the truncated 
system, and thus a local center-stable manifold for (KG) at time t = 0. 



2. We construct a global center manifold W c for the truncated system by constructing the solutions 



that are close to the origin for all times. (Subsection 4.3) 

The crucial point consists in linking the two previous steps: we prove that the solutions of the 
truncated system constructed in step 1 actually tend to some solutions contained in the local 
center manifold as t — > +oo and thus remain small for all times t > 0. As the central part 
remains unconditionally small, we conclude that these solutions are actually solutions of the 
original system (without truncation). This gives us a global center stable manifold W cs which 
is parametrized by the initial value V s of the stable part that we add to h and by V c the initial 



value of the elliptic part, both being small. (Subsection 4.4) 



For the proof of Corollary |1.3[ it remains to prove that we can choose V s and V c in such a way 
that the solution hits the reversibility plane b n = for all n: namely the solution at t = belongs 
to the reversibility plane as soon as V s = 0, and V c is symmetric. The corresponding solutions are 
then automatically symmetric and thus describe heteroclinic connections between two symmetric 



solutions of the center manifold. (Subsection 4.5) 



4.1. Preliminary results. 

by 



(4.1) 



X 



To begin with, we write an equivalent formulation of (|KG|). Denote 




A 



A + TO 2 



F(X) 



-u 





2p+l 



The equation ( KG ) is equivalent to 

(4.2) d t X = AX + F(X). 

By the result of Theorem |1.1| the equation (4.2) admits, for a given initial condition X(0) € X a 
unique global solution X G X so that sup||A(t)|| < C(||X(0)||#). We fix Co > large enough, so 

that all the initial conditions X(0) we consider in the sequel satisfy ||X(0)|| < Cq. Therefore all the 
solutions X we are going to consider satisfy 

sup\\X(t)\\<C(C ). 
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61, ct n , b n , ...) (schematic) 




Figure 3. Schematic overview of the strategy and the notations. 

We denote T> c the projection on X c of such initial conditions: 

V c :={QX\ \\X\\ < Co}. 
Let 6 e C°°([0, 00), K) be a cut-off function so that 

1, s < 5 
0, s > 25, 



e(s) 



5 being a small parameter smaller than 1. 
For any function G : X — > R, denote by 



G(X) = G(X6{\\X C \ 



Lemma 4-1- — Consider the function F given by (4.1) and let X,X' G X. Assume moreover that 



\Xh\, \ X' h \ < 5 < 1. Then there exists a constant C > such that 

(4.3) \\F(X)\\<C5 3 , 

(4.4) \\F(X)-F(X')\\<C5 2 \\X -X'\ 



Proof. — Writing X 



and the Sobolev embedding 



and u = a + U with a £ X^ and U £ X c we have by the fact that \a\ < 5 



\\F(X)\\ < ||(l«l + l^l^(ll^kO) 2P+1 || i2( M) 

< cs 2 p +l + c\\u\\fil +1) e(\\u\U) 

< C5 2p+1 + C\\U\\ 2 ^ 1 0(\\U\\ H i). 
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This gives (O), since \\U\\ H i6(\\U\\ H i) < CS. 
To get (4.4), just remark that for \u\, \u'\ < C5 



F ( U \ _ F H \ || 2 < J | n 2 P+ l _ n '2 P+ l|2 dx < J^ u _ u>l 2 d ^ 



and then apply this inequality to 



I =X0(||X c ||)and l 



X'0(\\X'\ 



□ 



We write the linearisation of (4.2) around the homoclinic orbit h = 
define 

(4.5) M{Z) := F{Z + h) - F(h) -dF[h].Z, 
with 

dF[k]= ( -(2p+°i)a 2 f(t) 
As a result, Z satisfies the equation 

(4.6) 3 t Z = C\Z)+M{Z). 



Set X = Z + h and 



Similarly as in Lemma 4.1 we have 



Lemma — Consider the function M defined in (4.5) and let Z,Z' £ X with \Zy\,\Z'A < 5 and 
with \\Z + h\\,\\Z' + h\\ < Ci 

\\M{Z)\\<C5\ 



(4.7) 
(4.8) 
(4.9) 



\\K(z)\\<c\\zf, 

\\N{Z)-K{Z')\\ <C5\\Z-Z'\\. 



4.2. The local center-stable manifold. — The equation (|4.6| is equivalent to the system 
(4.10) 



d t z h = c{z h + PM(z), 

d t Z c = ClZ c + QN(Z). 



We will be interested in solutions to (4.10) for which the center part is small, that's why we introduce 
the following truncated system 

f d t z h = c t h z h + PN(z), 
1 d t z c = CiZc + QNiZ). 



(4.11) 



We now look for solutions to ( |4.11[ ) with small hyperbolic part. Fix < r < m and denote by 

£+ = {Z £ C(R+,X), \\Z\\ r := supe- rt ||Z(t)|| < oo}. 

We introduce a new small parameter e which will quantify the size of Z, namely \\Z\\ < Ce 2 . We will 
choose e = <5 2//3 in such way the truncation 9 can be removed. 
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Proposition 4-3. — Let V c G X c with \\V C \\ < e 2 and let V s G E urai/i |T^| < e 2 . T/ien there exists a 
unique solution Z = Zy c y a to (4.11) in 

B+ = {Z££+, sup \Z h (t)\ < 5}, 
t>o 

and so that 

(4.12) <Z ft (0),<r*(0)) = V s , Z c {0) = V c . 
Moreover, this solution satisfies 

(4.13) sup|Z h (i)| < Ce 2 , 

t>o 

and 

(4.14) \\Z C \\ < Ce 2 , for all 0<t< e~ l . 

This result means, that for any fixed [|V^|| < e 2 , \V S \ < e 2 , there exists a unique choice of V u = 
{Zh(0), p*(0)) G K so that the corresponding solution has a small hyperbolic part for all times and 
a central part small for relatively long time (comparing with the time necessary in order that the 
homoclinic orbit h reaches a neighbourhood of radius e 2 of the origin). Let us call local center-stable 
manifold the set of all the solutions verifying this property: 

Definition 4-4- — For ||T^,|| < e 2 and \V S \ < e 2 , denote by Zy c y a the function given by Proposition 
|4J3]and 

x v c y a '■= z v c y 3 + h. 

We define W cs by 

W cs = |J {X Vc , Vs (0)}- 

\\V c \\<e 2 ,\V a \<e 2 

Theorem 4.10 below will ensure that for e sufficiently small, this set is actually the standard center- 
stable manifold of the equilibrium for the initial equation (4.2). 

Proof. — We define the map T : — > 

TZ(t) = V s a(t)+K(t,0)V c + f (PAL(Z){r),a*(T))dTo-(t) 

Jo 

30 rt 

{PM(Z)(T),p*(r))dT P (t)+ / K(t,T)QK(Z)(T)dr, 



and show that for 5 > small enough, J 7 is a contraction in Bf. 

• To begin with, we show that T : Bf — > Bf . Since P = J M , we have 

i<PAa^*>i < \pm(z)\w*\ < \\mz)\\L*(M)e mT < \\mz)\\e % 

Let \Z h \ < 5, then by ( |47| , for all t > 

"t /"+oo 



/•t r+oo 

\(P(Z)) h (t)\ < C\V s \e- mt + C \\N{Z)\\e mT &Te- mt + C / ||Ar(Z)||e- mT dre 
(4.15) < C(e 2 + <5 2 ) < Ce 2 . 
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On the other hand by Proposition 3.4 and (4.7) 
(4.16) 



\\(F(Z)) c (t)\\<C\\V c \\ + C [ \\AL(Z)\\dT<C(\\V c \\+5 2 t), 

Jo 

which in turn implies that ||(j r (Z)) || r < oo. 



We now show that J 7 is a contraction in Bt. By Lemma 
\{J : (Z 1 )) h (t)-(T(Z 2 )) h (t)\< 



4.2 



< C / \\K(Zi) - K(Z 2 )\\e mr dTe- mt + C / \\M{Zi) - M{Z 2 )\\e- mT dre mt 

Jo Jt 

ft f + OO 

< C8e- mt \\Z x -Z 2 \\ r / e( m+r ) r dr + C8e mt \\Z x - Z 2 \\ r / e(~ m+r ) T dr 

Jo Jt 



and therefore 

(4.17) 

Similarly 



and then 
(4.18) 



\{HZi)) h {t) - {T(Z 2 )) h (t)\\ r < CS\\Z X - Z 2 \\ r . 



\{F{Zi)) c (t) - {HZ2)) c {t)\\ < C mZi) ~ K(Z2)\\dr 

Jo 

< CS\\Z X - Z 2 \\ r [ e rr dr, 
Jo 



\{HZi)) c (t) - (nz 2 )) c (t)\\ r < C8\\Z X - Z 2 \\ r . 



Thus, ( |4.17 ) and (4.18) show that J 7 is a contraction whenever < 5 < 1 is small enough, and we can 
deduce that there exists a unique fixed point Z G . By definition of J 7 , it is clear that this solution 



satisfies (4.12) 



Finally, in view of the choice e = 5 2 ^ 3 , the bound (4.13) comes from (4.15) and (4.14) is a direct 
consequence of (4.16). □ 



4.3. The local center manifold. — In this subsection we want to construct solutions that remain 
close to the origin for all time. Actually this will be achieved by constructing a local center manifold. 



Again we consider the approximation of the initial problem (4.2 ), in which we truncate the non-linearity 

f d t X h = A h X h + PF(X), 
(4 ' 19) I d t X c = A C X C + QF(X). 



The next result shows that there exist solutions to (4.19) with small hyperbolic components. For 
< r < m we define 



£ r = {X € C(R+,X), \\X\\ r := supe- r l*l||X(t)|| < 00}. 
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Proposition 4-5. — Let V c 6 T> c . Then there exists a unique solution to (4.19) in 

Bs = {Xe£ r , suj>\X h (t)\<5}, 



and so that 
(4.20) 

Moreover, this solution satisfies 
(4.21) 



QX(0) = V c . 



sup\X h (t)\ < C5 3 . 

teR 



The set of these solutions is the center manifold of for equation (4.19) : 



Definition 4-6. — For V c £ T> c , denote by Xy c the function given by Proposition 4.5 We define W c , 
a global center manifold for (4.19) (and a local center manifold for (KG}), by 



w °= U {XvMY 



V c eV c ,teR 



P0 



4.3 



Denote by oq 



1 



??? 



Proof. — We proceed similarly as in the proof of Proposition 

^ \ . The dual basis of {a , po} reads {<rg, p$} with = - ^ I 1 J ana p§ = 2 , 

It is clear that (doe - " 1 *, poe mt ) form a basis of the solutions of dtX^ = A^Xh- Denote by 
propagator of the equation dfX c = A C X C . Then for X £ E r we define 



and /Jq — 



and 

ni 
1 

K the 



(4.22) 



QX{t)=K{t)V c + I (PF(X)(T),a* e mT )dTo- e- 

J — oo 



■mt 



r+oo ft 

/ (PZ(X)(r),^e-^)dr /O0 e m *+ / K{t - r)QF{X){r)dr. 
Jt Jo 



It is straightforward to check that any fixed point of Q in B$ satisfies (4.19) and QX(0) = V c . 



Firstly we show that Q : B$ — > B$. Let X £ B$, then by (4.3) 



/t f+OO 
\\F(X)\\e mT dTe- mt + C / \\F(X) || e - mr dre m * 
-OO Jt 



(4.23) 



Next, by (4.3) again 



\{Q{X)) c (t)\\<C\\V c \\+C f \\F(X)\\d T < C + C8 3 t, 

Jo 
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which in turn implies that ||(£7(JT)) || r < oo. Hence, with (4.23) we get Q{X) G Bg. 
• In a second time, by (|4.4|) we can write 



\(g(x 1 )) h (t)-(g(x 2 )) h (t)\< 

/t r+oo 
WKX,) - F(X 2 )\\e mr dTe- mt + C / - F{X 2 )\\e- mT dre mt 

-oo J t 

/t /*+oo 
e mr+r\r\ dT + C 8 2 e mt \\Xi - X 2 \\ r / e- mT+r l r ldT, 
-oo J t 



and therefore 

(4.24) 

Similarly 



" {g(X 2 )) h (t)\\r < CS 2 ]^ - X 2 \\ r . 



\{G(X 1 )) c (t)-(g(X 2 )) c (t)\\ < C \\F(X 1 )-F(X 2 )\\dr 



< C<5 2 ||Xi -X 2 \\ r / e r l T ld7 



and then 
(4.25) 



|(£(Xi)) c (i) - {g(x 2 )) c (t)\\ r < cs^x, - x 2 \\ r . 



Thus, (4.24) and (4.25) show that g is a contraction in Bg whenever < 5 < 1 is small enough. The 
bound tell is given by 14.231). □ 



We are now able to give a parametrisation of the local center manifold given in Definition 4.6 Let 
\& : D c — > Xh the map defined by 

*(V C ) = (X Vc ) h (0), 

then 

(4.26) W c = {(*(y c ),V e ), V c G V c }, 

and, in view of (4.22), ^ is quadratic at the origin, ^{V c ) = 0{\\V C \\ 2 )- This is a particular case of a 
result of Mielke [21^ 

As a consequence we can prove that the origin is Lyapunov stable within the center manifold: 



Lemma 4-7. — Let X be a solution of (4.19) which lies in W c . Assume that for some time t* > 0, 
\\X c (t*)\\ < Ce 2 . Then for all time t G R 

(4.27) \\X c {t)\\<Ce 2 . 



Proof. — It is here convenient to work in the coordinates (a n ,b n )n>i- Recall that from (2.6) the 
Hamiltonian of (4.2) reads 

H° = i(6g - m 2 a 2 ) + J + 2 ^ £ (a + U) 2p+2 dx, 
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and obser ve th at J/C < \\X C \\ 2 < CJ. We first show that \H°\ < Ce 4 . By assumption, J(t*) < Ce 4 , 
while by (4.21) we obtain 

(4.28) \ao(t)l < C6 6 = Ce 9 < Ce 2 for all t G M. 

Next, using Lemma 2.2 we obtain, for t = t* 

which proves the claim. Next, with (4.28) we can write for all t G R 

J(t)<H°-~(b 2 -m 2 a 2 )<Ce\ 



which implies (4.27). 



□ 



4.4. The local center-stable manifold is global. — Let \\V C \\ < e 2 , \V S \ < £ 2 and consider the 
solution Z = Zy y a to (4.11) given by Proposition 4.3 Set 

(4.29) 



Ami. 

m e 



By (1.14), with this choice we have \h(t)\ < Ce 2 for all t > t £ /2. Next, we define 6 G C°°(R,R) so that 



0, t < t £ /2 

1, t > t F , 



and with \0'(t)\ < C/t e for all t G R. We set 

(4.30) Y(t) = Q(t)(Z(t) + h(t)). 



Proposition 4-8- — Consider the function Y defined in (4.30) . Then there exists X solution of (4.19 ) 
within W c so that for alltER 

(4.31) \\(Y-X)(t)\\ < Ce' 2 e- rt and \{Y - X) h {t)\ < Ce 2 . 
Proof. — The function Y satisfies the equation 

d t Y = AY + F(Y) + 5, 

with 

(4.32) S(t) := e'(t)(Z + h){t) + @(t)F((Z + h)(t)) - F{0(t)(Z + h)(t)). 

Consider now a solution X of the equation dtX = AX + F_{X) and set A = Y — X . Then A satisfies 

(4.33) d t A = AA + F(Y) - F(Y - A) + S. 



The proof of the proposition consists in finding a solution of (4.33) so that (4.31) holds and so that 
X(0) := y(0) - A(0) G W c . 
Set 



£~ = {A G C(K,#), ||A|| r := sup e r *||A(i)|| < oo}. 



and 



{Ag«S t 7, sup I A h (t) I <5}, 
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where again we choose S 2 = e 3 . 

We claim that we can define a contraction J : Bj — > Bj by 

JA(t) = f (P{F(Y) - F(Y - A) +S)(T),a* e m -)draoe- mt 

J — oo 

/+oo 
( P ( FiX) -F(Y-A)+S) (r) , p*e-™ )dr Po e mt 

/+oo 
K(t - t)Q(F(Y) - F(Y - A) +5)(r)dr. 

It is clear that the fixed point will satisfy ( 4,33[ ). In a second time we will check that (4.31) holds. 
• We begin by estimating 5. First observe that by Proposition 4.3 

(4.34) \\Z(t)\\<Ce 2 for all < t < e~ l . 

Then, for e > small enough we have t £ < e^ 1 and we can write for all t 6 [t e /2, t e ] 

(4.35) \\Z(t)\\ < sup \\Z(t)\\e- r{ t-^ < CeV^^ < Ce~ 2 e 

te[t s /2,t e ] 



-rt 



since e r e < e m e = e . Denote by \ the indicator of the interval [t £ /2,t £ ]. As Z + h is small on 
[t e /2,t e ] we deduce from (4.32) that 

(4.36) 



C 



S{t)\\ < r (\\Z(t)\\ + \\h(t)\\) X (t) + C(\\Z(t)\\ + \\h(t)\\) X (t). 



As a result, from (4.35) and (1.14) we infer 

(4.37) \\S(t)\\ < C(e- 2 e- rt + e- mt ) X (t) < Ce-\- rt X {t). 

The inequality (4.36) together with ( 4.34[ ) also gives 

(4-38) mt )\\<C-^- X {t). 



• We have |Y/j| < 5. For A G B s we can apply (4.9) to deduce 

\\(F(Y)-F(Y- A))(r)|| < C<5 2 ||A(r)|| < C6 2 \\A\\ r e- rT . 
Thanks to this latter inequality and ( 4.37[ ), we obtain 

llJA^He^ < C f (<5 2 ||A|| r + e- 2 x(r))e( m - r ) T dre('- m ) i 

J — oo 

-+oo 



-.rt 



and therefore 
(4.39) 



/+0O 
(5 2 ||A|| r + £- 2 x(r))e-( m+r ) T dre( m+r ')* 

/+oo 
(5 2 \\A\\ r + e- 2 X (r))e- rT dTi 

< C(5 2 \\A\\ r + e- 2 ), 

\\JA\\ r < C{5 2 \\A\\ r + e- 2 ). 
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Next, by fl4.38| ), for all t G R 

|(^A) h (t)| < C/* ( < 5 3 + T -^ x (r))e-dre 

J —CO I ^ ^ I 

+C Jf (5 3 + _ x ( T ))e-^ dr 
(4.40) < C(<5 3 + e 2 ). 



mi 



The bounds (|4.39|) and (I4.40I) show that J : 07 



• Let Ai, A 2 G 07. By ( |4.9| ), we have 
(4.41) 
Clearly, 



\\(F(Y - A 2 ) - F(Y - A!))(r)|| < C^Ai - A 2 || r e 
JAi(t)- JA 2 (t) 



/ ( P( F(Y - A 2 ) - F(Y - Ai), a*e mT )dr ae" mt 

J — oo 

/+oo 
( P( F(y - A 2 ) - F(Y - At) ) (t), p*e~ mT )dr pe 

/+oo 
if (i - r)Q( F(Y - A 2 ) - F(Y - A x ) ) (r)dr, 



and by fl4T4Tj ), for all i G R 

||JAi(t) -JA 2 



< C f (5 2 ||Ai - A 2 || r )e (m - r)T dre( r - m )* 

«/ — oo 

/+oo 
(<5 2 ||Ai - A 2 || r )e-( m+r ) r dre( m+r )* 

/+oo 
(<5 2 ||A 1 -A 2 || r )e- rr dre r 



< C5 2 ||Ai - A 2 || r . 

As a consequence, for 5 > small enough, the map J7 : £>7 - 
unique fixed point A G 07. By (4.40) and choosing again 5 2 = e 3 , we get |Afc(i)| < e 2 . Furthermore, 
as A = J (A), (|4.39|) leads to ||A|| r < Ce~ 2 , which in turn implies ||A(t)|| < Cs~ 2 e~ rt for all t G R. 



5 is a contraction, thus there exists a 



• We now define X = Y — A and it remains to show that X G W c . By definition, it is sufficient to 
prove that \X h (t)\ < S for all t. Write Y = Q(Z + h). By ( |4.13| ) and the choice of t E we get that 
lift (*) I < Ce 2 for all t G R and then 

\X h (t)\ <\Y h (t)\ + \A h (t)\ <Ce 2 <5, for all t G R, 

which implies the result, for e > small enough. □ 



Proposition 4-9- — Let \ V s \, \\ V c \\ < e ande small enough. The solution Z of (4.11 ) with Z c (0) = V c 
and Z s (0) = V s constructed in Proposition 4-3 satisfies \\Z c (t)\\ < 5 for all t > 0. In particular 
X = Z + h is a solution of (4.2 ) for t > 0. 
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Proof. — (see Figure [3] on page 
Proposition 4.8 For t > t e , we have Y = Z + h, thus by (4.31 ) 



15) For such a Z, consider Y defined by (4.30) and X defined in 



\\X c (t) - Z c (t)\\ = \\X c (t) - Y c (t)\\ < \\X(t) - Y(t)\\ < Ce~ z e- rt , for t > t e . 
4 1 

Recall that by (4.29), t £ := — In -. Thus, for e > small enough t £ < e -1 and the previous inequality 
1 ' me 

implies 



(4.42) 



\X c (t) - Z c (t)\\ < Ce~ 2 e- r£ 1 < Ce 2 , for t > e~ x . 



Now, take t* = e~ x , by Jili]) we have ||Z c (t*)|| < Ce 2 , and then ((442]) implies ||X c (t*)|| < Ce 2 . By 
Lemma 



4.7 



we deduce that ||X C (£)|| < Ce 2 for all t > e 1 , and coming back to (4.42), we infer 

< Ce 2 < S, for all t > e~ X . 



This bound together with (4.14) conclude the proof. 



□ 



Gathering the results of the previous propositions, we get the following Theorem, which ensures that 
W cs is the standard center-stable manifold of for equation (4.2). 

Theorem 4-10. — For e sufficiently small, for all <[ \V C <[ \,\V S \ < e 2 , the functions 

Xv c ,Vs = z v c ,v s + h, 



(where Zy c y s is given by Proposition 4-3) defined for t > are solutions of (4.2). Moreover, Zy C) y a 
belongs to and satisfy 

\\Z c (t)\\ < Ce 2 < 5, \\Z h (t)\\ < Ce 2 , for all t G M + ; 



and there exists X G W c such that 

\\(X Vc , Vs -X)(t)\\<Ce- 2 e- rt . 
This completes the proof of Theorem |1.2| 



4.5. Construction of reversible heteroclinic solutions: proof of Corollary 1.3 — Recall 
that the vector field is reversible, and thus if the initial condition of a solution X satisfies 

X(0) = SX(0), 

then the solution is reversible, i.e. 

X(-t) = SX(t), for all t G M. 
Lemma 4-H- — Fore sufficiently small, for all \\V C \\ < e 2 verifying 

v c = sv c , 



Xy c y s defined in Theorem 4-10 satisfies 



X Vc , (0) = SX Vcfi (0). 
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Proof. — On the one hand, from the explicit form (1.13) of h, we know that h(0) = Sh(0). On the 
other hand, given that Z is a fixed point of the map T defined in the proof of Proposition |4.3[ we get 



/■+00 

Z Vci0 (0) = V c - / (PAA(Z)(r),/(r))drp(0). 
Jo 



And since p(0) 
SZ VctO (0). 



(see Lemma 



3.5), we have Sp{0) = p(0). So, if V c = SV C , then Z Vcfi (0) 



□ 



As a consequence of the previous lemma, Xy c fi is a reversible solution. From Theorem 4.10, we then 
get that on the one hand, Zv Cl o = Xv C7 o ~~ h belongs to £ r and satisfy 

\\(Z Vc ,o) c (t)\\ < 5, \\(Z Vc , ) h (t)\\ < Ce 2 , for all t G R; 

and on the other hand, there exists X G W c such that 



\\(X Vcfi -X)(t)\\<Ce- z e 



2„-rt 



for all t > 0, 



and 



||(AV cj0 - SX)(t)\\ < Ce~ 2 e rt for all t < 0. 



This means that Xy ct o is an heteroclinic connection between the two solutions X and SX of the center 
manifold W c . 

Then Corollary 1.3 follows since the condition V c = SV C and V s = is equivalent to (1.17). 



5. Solutions close to the homoclinic orbit 



When we impose a n — b n — for xi ^ 1 equation (|KG|) reduces to (11.121) and this two-dimensional 
system admits periodic orbits (denoted by JC V in Section 1.5) "inside" the homoclinic orbit h (see 



Figure [TJ. In this section we prove Theorem 1.5 which asserts that, when turning on the modes n > 1 
(but with a small amplitude) we still have solutions close to the periodic orbits of (1.12). 



1.5 



i . 



5.1. Proof of Theorem 

sition u = aoeo + U, the Hamiltonian H in (jl . 1 0|) reads 
(5.1) 



Recall the notations (2.2) and (2.3). Then with the decompo 



H=l(bl-rn^) + J +1 ^J M (a Q + U) 



2p+2 



dx. 



We choose initial conditions so that ao(0) = T], & (0) = and J(0) < rf with < rj < 1. Then 



(5.2) 



H 



m 2 r] 2 + 0{rf 



Lemma 5.1. — For rj small enough, as long as J < Ci] 7 / 2 , there exists C\ = C\(m) > such that 
(5.3) 77/2 < a < d, |6 | < C x . 



26 



BENOIT GREBERT, TIPHAINE jEZEQUEL AND LAURENT THOMANN 



Proof. — Assume that J < Crj 7 ^ 2 . First we prove that ao cannot vanish: Indeed if it was the case, 
by (5.1) we would have H > which is in contradiction with (5.2). Next we prove that there exists 



C > so that 
(5.4) 



2(P 



Write the binomial expansion 



JM ~[ V 3 J J A 



(5.5) 



M 



Apply Lemma 2.2 and use (2.5 ) to each term of the sum with e > small enough, and get 
2p+2 , N 2p+2 

V( P Vo P+2_j / U*dx <C Ya 2 p+2 -iJ^ 2 < \ a 2p+2 + CJ p+ \ 
~[ V 3 J Jm \ 1 



which together with (5.5) yields (5.4). Now, for r\ > small enough, CJ P+1 < J and by (5.1) and 

1 



(5.4) we obtain 



(5.6) 



1 



2^~"' ^ 4(p + 1) "° 

from which we deduce that |6q| < C\ and < ao < C\ for some constant C\ depending on m. 
Furthermore, the bound (5.6) implies —m 2 a^/2 < H < — |m 2 ?7 2 which completes the proof 

Lemma 5.2. — Assume that J < Cifl 2 . Then there exist times < T v < T v so that 

|ao(T„) - (p + if^m 1 '^ < Cri 7 ' 2 , b (f v ) = 0, 

\ao(T n ) - ri\ < Crf/ 2 , b (T v ) = 0, 
and so that bo > on (0,T„) anc? 6o < on (r„,T„). 

Proof. — Recall that (ao,&o) satisfy the equations (the (afc)/j>i are here considered as source terms) 



m a Q ) + 



1 



2p+2 



< H 



1 



-m 2 rf + C(r/ 3 ) < -^m 2 r] 2 , 



□ 



(5.7) 



d = b , b = m a 



[ (a + U) 2p+1 dx. 
Jm 



Since )/ < 1, we have &o(0) > 0, and thus for small times t > 0, ao(t) > and bo(t) > 0. We now show 
that there exists > so that 6p( T r? ) = and bo(t) > on (0,^). By contradiction, if &o > for 
all t > 0, from the first equation in (5.7), we get that ao is increasing, and thus, since ao is bounded, 
ao(t) — > L > when t — > +oo. 

Next, we expand the second equation in (5.7) (using that / Udx = 0) 

Jm 

b = m 2 a -a 2 p+1 + O(r ] 7 / 2 ), 



and therefore when t 
(5.8) 



+oo 



b = m 2 L-L 2 P +1 + 0(ij 1 ' 2 ). 
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Since bo is bounded, there can not exist any constant cq > so that &o > c or ^0 < ~ c when 
t — > +00. Thus from ( |5.8[ ) we deduce m 2 L — L 2p+1 = 0(rj 7 ^ 2 ). The implicit function theorem then 
gives that L = 0{r] 7 / 2 ) or L = m l / p + 0(r] 7 / 2 ). We can exclude the first case since L > ao(0) = 77. We 
now use the conservation of H: from (|5.1l) and the fact that bo > we get when t — > +00 



bo(t) 



P+l 



-m 



This is in contradiction with do = 60 and the fact that ao is bounded. As a conclusion there exists a 
smallest time T v > so that 60 (T„) = 0. 

With the conservation of Jf we find a (T v ) = (p + I) 1 /! 2 ? 3 )™ 1 /? + C>(r/ 7 / 2 ). 

Similarly, we show that there exists > so that 60 < on (T„,T„), bo(T„) = and ao(T^) = 
?? + 0(r? 5 / 2 ). ' ' □ 



Denote by 



I := 



^ 00 



2 2 
fc -m 



+ (2p + l)a^)4 + 6i 



fc=i 



(5.9) 



I < 



Lemma 5.3. — Le/j < T v < T v be given by Lemma 5.2 and assume that J < Crfl 2 . Then we have 
the a priori bound 

Ca^aal + CI 3 / 2 if t e (0, f n ), 

. CI 3 / 2 if te(T v ,T v ). 

Proof. — We write u = aoeo + U, then for k > 1 

/ u 2p+1 (-,x)e k (x)dx = (2p+ l)a 2 p a k +Pk, 
Jm 

where pk is defined by 

(5.10) p k = j ((a eo + U) 2p+1 - a 2p+l - (2p + l)afu) e k dx. 



Thus the equation on (a k , b k ) reads 
(5.11) a k = b k , b k 



m 



+ (2p + l)a 2p )a k -p k . 



Observe that since ao is bounded, there exists C > so that for all (a k , b k ) k >\ G h 1 x 
(5.12) J<I<CJ. 



Now with (5.11) we compute 
(5.13) 



/ = p(2p + l)ap p ^0 ^2 a\ - ^2 b kPk < p(2p + 1)oq P ^0 ^ af, + 

k=l k=l k=l 



2 (N*)\\P\\£ 2 (n*)- 
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We have 



\i 2 CH*) 



< CI 1 / 2 . Then by Parseval and the fact that |ao| < C 



2 (N*) 



\(a + U) 2p+1 



2p+l 



(2p + l)a 2 p U\\ L2 



(M) 



< C||LT 2 || L2(M) + C||^ +1 || i2(M) 



c 'II c/ IIl 4 (m) 



T^\\ U ll i 2(2 P +l)( M ) 



< CI + C/ p+ 2 < CT, 

Denote by 6 (0,1^) the time so that bo{T fj ) = 0. Then 



2.2 



(5.14) 

where in the last line we used Lemma 
since do > on (0,T^) and do < on (1^,7^), by (5.13) and (5.14) we deduce (5.9). 

Lemma 5-4- — Assume that J < Crj 7 ^ 2 , then the following bounds hold 



□ 



(5.15) 



1 



cm- <T„< Cin 



1 



: hi - < T„ < Cin-. 
r\ r) i] r\ 

Proof. — • We focus on the control of T v , the control of is obtained similarly. Then fix 5 := m 1 / p /2 
and let < t v < T v be so that a (T v ) = 5. Then since \J clq + o^ > c on [t v ,T v ], we g et Trq T<q ^ C 
Hence it is sufficient to prove that c In - < t„ < C In - . 
• Denote by 

(5.16) /(x) = -mV + 

Let ip be the function defined by 

i?l\{t) = -2(J(t) - J(0)) 



1 



p + 1 
1 



-X 



2p+2 



P+l 



(po(i)-Po(O)). 



By (5.1) and (5.10), tp satisfies 

(5.17) &§ + /(ao) = /(»j)+»7 7/ V(<)- 
Let us prove that, as long as J < Crfl 2 and t < r„, 

(5.18) < Cmin(l,i 2 ). 

Since J < Crfl 2 and by (5.12) together with ( |5.14 ), we have \ip\ < C. Then, since &o(0) = and 
do(0) = &o(0) = 0, from ( |5.17 ) we infer that 0(0) = 0. To obtain (5.18), it is therefore enough to prove 
that \(p\ < C. This is done by computing J and po with the relations (5.11 ), and with Cauchy-Schwarz 
we can check that | J| < Crj 7 / 2 and |po| < Tp 7 / 2 . 
• Since on [0,t^], bo > 0, from (5.17) we get 

ao = b = ^f(v)-f(ao) + r] 7 /Mt), 
and by integration and the change of variables a = ao(t) we obtain 



ao(t)dt 



da 



o Vf(v)~ f(a ) + v 7/2 v(t) h ^//(?7)-/(a) + 7? 7 /2^( a -i( a )) 
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• We now show that if 5 > is small enough (but independent of rj) there exist c, C > so that for all 
r] < a < 5 



(5.19) 



V ~ rj 2 ) < f( V ) - f(a) + r, 7 l\{t) < C{a 2 - if), 



and this will complete the proof, since 



da 



\/ a 2 — if 



In 



It is clear that c(a 2 — rj 2 ) < f{rj) — f(a) < C(a 2 — rj 2 ), it is then enough to show that rf/ 2 ip is a 
negligible amount when rj -C 1. From (5.7) we deduce that for t G [0,tJ, clq > ci] and by Taylor, 
\ao(t) — cio (0)| > crjt 2 which in turn implies 

C 

\% 1 ( a )\ = \ a o V) ~ a o 1 ( r l)\ < — I « - V\ 1/2 - 



Then by (5.18), for r\ < a < 5 

7? 7/2 |<^(ao \a))\ < Crfl 2 rr 2 \a - rj\ < Crf' 2 {a 2 - rj 1 ), 



which proves (5.19). 



□ 



We are now able to complete the proof of Theorem 1.5 (i). In view of Lemma 5.2 it suffices to prove 
that I < Crj 1 / 2 for t 6 [0, T„] (recall that J < I). This is achieved by a bootstrap argument: In the 
regime t € (0, T„), as long as / < Crfl 2 , we have I < C(aoa 2 ] p ~ 1 + rf^)I which implies 



Krf 



lit) < I e c fo v (*o(s)% p 1 (s)W /A )ds < KIq 

for some absolute constant K > 1 since a® < C and rf'^T„ < 1 for r\ small. In particular I(T„) < Krf. 
Next, in the regime t £ (T„,T„), and as long as 7 < Crfl 2 we have 7 < Crf^I which implies by 
integration for all t G (7^, 7^) 

(5.20) I(t) < /(f^e^ 774 ^) < 2KI(f v ) < 2Krf . 

Using the conservation of the energy again, we then get that \ao(T v ) — (p + l) 1 /( 2 P)m 1 / p | < Crf as 
well as \ao(T ri ) — 7]\ < Crf, (which is slightly better than the result of Lemma 5.2). 



5.2. Proof of Theorem 1.5 (ii). — Consider a trajectory of (1.11) with ao(0) = rj, bo(0) = and 
J(0) < rf. Thank to (5.20) we see that we can iterate the argument as long as 7 ^ Crj' 1 1 2 ', namely 



In-) times (each time corresponds to a loop). This gives that for \t\ < (In 1 ) , J < Crf I 2 . 



l\2-<5 



Next, consider the trajectory of (1.12) described by (a^,/^) 
the definition (5.16) of /, then the conservation law for (a^, (3 V ) 

% + /(«„) = f( V ). 
5 there exists ^a v ,^„) 



so that a v (0) 
reads 



rj, p v (0) = 0. Recall 



(5.21) 

We now claim that for each \t\ < (in 



which satisfy (5.21) so that 



\OLr 



a (t)\<Cr, 5 / 2 , \(3 v -b (t)\<Cri 5 / 2 
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• Fix X p G (l, (p + l) 1 '^). Assume that ao ^ [X p l m l /' p , A p m 1/,p ]. Then we define (3^ = bo and from 
( Pl| ) 



Next from (5.17) we deduce 



and since |/'(ao)| > cqT] we can apply the mean value theorem to get 

\a v -a \ < ^l^lsupK/- 1 )'! < Cr? 5 / 2 . 

• In the region ao £ [A~ 1 m 1 / p , A p m 1//p ], we define g(bo) = b^/2. Then we have |</(&o)| > c i an d we can 
perform a similar argument to get a v = ao and — 6q| < C?? 3 . 



Appendix A 

In this appendix we recall a result concerning the long time behaviour of the solution of the ordinary 
differential equation 

x = (-a 2 +q(t))x, f G 1 

where a is a real constant and t \— > q(t) is a continuous function. Of course the solution of such 
linear equation are globally defined, but we would like to know whether the solution are bounded or 
not. It tur ns out that even if ||o||oo < ° 2 an d q(t) — > when t — > ±oo, the solutions may grow 



indefinitely ^ The good condition concerns the integrability of \q\: 



Lemma A.l. — Assume that / \q(t)\dt < +oo then all the solutions of the Cauchy problem 

Jr 



(A.l) 



x = y 

,2 



[y = (-a +q(t))x, 

are bounded on EL More precisely, there exists C > which only depends on q so that 

(a 2 |x(*)| 2 + \v(t)\ 2 ) < C(a 2 \x(0)\ 2 + \y(0)\ 2 ) for all t G R. 

The proof is classical (see for instance [14, page 212]) but does not precise the bound on the solution 
in term of the initial datum. The following argument is more explicit. 



Proof. — Denote by S the flow of ( |A.l | for q = 0, and introduce the norm ||x,y|| = (a 2 x 2 + y 2 ) 1 ^ 2 



(observe that S preserves this norm). For T > 0, the Duhamel formula reads 

(A - 2) ( m ) = s(t - T) ( m ) + SI s «- T -4 gW ° (,) l d »- 



2. Consider the example induced by the solution x(t) = sin£(l + (t — s'mt cost) 2 ) (see O. Perron |23| ). 
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f Tl 

There exists T\ so that / |g(i)|di < 1/2. Therefore from (A. 2) with T = we deduce that for all 

Jo 

t G [0, T\] 

\\x(t),y(t)\\ < ||x(0),y(0)|| + \ max \\x(s),y(s)\\, 

2 sG[0,Ti] 

and thus max ||x(s),y(s)|| < 2||x(0),y(0)||. 

se[o,Ti] 

By induction, we can define a finite number of times T2, . . . , so that 

"+00 



-j+i 



\q(t)\dt = 1/2 and 



r+00 

/ |g(i)|df < 1/2. Then we apply (A. 2) with T = Tj and we show max ||x(s), y(s)|| < 
•//, " se[Tj,Tj +1 ] 

2\\x(Tj),y(Tj)\\ which in turn implies 

\\x(t),y(t)\\ < 2 k+1 \\x(0),y(0)\\ for all t > 0. 

□ 
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